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The  presence  of  Anopheles  (Nyssorhynchus)  dunhami  Causey  in  Colombia  (Department  of  Amazonas)  is  con¬ 
firmed  for  the  first  time  through  direct  comparison  of  mtDNA  cytochrome  c  oxidase  I  (COI)  barcodes  and  nuclear 
rDNA  second  internal  transcribed  spacer  (ITS2)  sequences  with  topotypic  specimens  of  An.  dunhami  from  Tefe, 
Brazil.  An.  dunhami  was  identified  through  retrospective  correlation  of DNA  sequences  following  misidentification 
as  Anopheles  nuneztovari  s.l.  using  available  morphological  keys  for  Colombian  mosquitoes.  That  An.  dunhami  oc¬ 
curs  in  Colombia  and  also  possibly  throughout  the  Amazon  Basin,  is  of  importance  to  vector  control  programs,  as 
this  non-vector  species  is  morphologically  similar  to  known  malaria  vectors  including  An.  nuneztovari,  Anopheles 
oswaldoi  and  Anopheles  trinkae.  Species  identification  of  An.  dunhami  and  differentiation  from  these  closely  re¬ 
lated  species  are  highly  robust  using  either  DNA  ITS2  sequences  or  COI  DNA  barcode.  DNA  methods  are  advocated 
for  future  differentiation  of  these  often  sympatric  taxa  in  South  America. 
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Anopheles  ( Nyssorhynchus )  dunhami  Causey  was  first 
described  from  Tefe,  Amazonas  (AM),  Brazil  (Causey 
1945),  but  was  placed  in  synonymy  with  the  malaria  vec¬ 
tor  Anopheles  nuneztovari  Gabaldon  (Lane  1953),  where 
it  remained  for  some  decades  until  taxonomic  revisions 
resulted  in  its  re-elevation  (Peyton  1993).  In  the  same 
study,  Peyton  also  reduced  the  widespread  east-Andean 
taxon  Anopheles  trinkae  Faran  to  a  junior  synonym  of 
An.  dunhami,  based  on  similarities  in  male  genitalia 
(Peyton  1993).  Later  integrated  studies,  combining  egg 
morphology,  cytogenetics  and  rDNA  sequences  of  the 
second  internal  transcribed  spacer  (ITS2),  confirmed 
that  An.  trinkae  and  An.  dunhami  were  indeed  separate 
taxa  and  An.  trinkae  was  reinstated  to  species  status 
(Lounibos  et  al.  1998).  Recently,  Calado  et  al.  (2008)  re¬ 
affirmed  the  differences  between  An.  dunhami  and  An. 
nuneztovari  using  the  morphology  of  male  genitalia  and 
DNA  sequences  of  the  ITS2  and  a  portion  of  the  mito¬ 
chondrial  cytochrome  c  oxidase  I  (COI).  Overlapping 
morphological  characters  between  proven  malaria  vec- 
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tors  An.  nuneztovari,  An.  trinkae  and  Anopheles  oswal¬ 
doi  (Hayes  et  al.  1987,  Olano  et  al.  2001,  Quinones  et  al. 
2006,  Gutierrez  et  al.  2009)  and  isomorphic  non-vector 
species  such  as  An.  dunhami  continue  to  hamper  correct 
species  identification  in  South  America  and  can  impact 
the  effectiveness  of  vector  control  programs. 

As  a  result  of  its  historically  confused  taxonomic 
status,  few  verified  distribution  records  exist  for  An. 
dunhami.  Records  to  date  are  limited  to  AM:  Tefe 
(type  locality)  (Causey  1945,  Peyton  1993),  Tabatinga 
(Lounibos  et  al.  1998)  and  Parintins  (Calado  et  al. 
2008).  Although  Tabatinga  is  adjacent  to  the  Colom¬ 
bian  border,  An.  dunhami  has  never  been  reported  from 
Colombia  and  the  species  is  not  included  in  the  avail¬ 
able  morphological  keys  to  Colombian  taxa  (Suarez  et 
al.  1988,  Gonzalez  &  Carrejo  2007).  Herein,  we  molec- 
ularly  characterize  topotypic  An.  dunhami  using  COI 
barcodes  [after  Hebert  et  al.  (2003)]  and  ITS2  rDNA 
sequences  and,  through  molecular  correlation  of  se¬ 
quence  data,  confirm  the  presence  of  An.  dunhami  in 
the  Colombian  Amazon  for  the  first  time. 

MATERIALS  AND  METHODS 

Mosquito  specimens  -  A  total  of  4,893  mosquitoes 
were  collected  in  four  departments  in  Colombia  (Antio¬ 
quia,  Amazonas,  Caqueta  and  Norte  de  Santander)  in 
2006  as  part  of  the  PhD  study  of  Ruiz  (2010).  Of  574 
mosquitoes  collected  in  Amazonas,  where  An.  nunezto¬ 
vari  s.l.  is  unreported,  two  specimens  were  keyed  out 
as  An.  nuneztovari  s.l.  using  the  morphological  keys  of 
Suarez  et  al.  (1988)  and  Faran  (1980).  DNA  was  extracted 
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TABLE 


Field  codes,  collection  localities  and  co-ordinates  of  specimens  of  Anopheles  dunhami  and  Anopheles  nuneztovari  C  used  in 
this  study,  alongside  GenBank  accessions  for  second  internal  transcribed  spacer  (ITS2)  and  cytochrome  c  oxidase  I  (COI) 
sequences  where  appropriate  outgroup  taxa  were  chosen  based  on  their  close  morphological  similarity  to  An.  dunhami 


Species 

Country,  locality 
(co-ordinates) 

Field  code 

GenBank 

ITS2  COI 

An.  dunhami 

Brazil,  Amazonas,  Tefe 

SA26 

HQ020394 

- 

(-03.3207,  -64.7235) 

SA44 

HQ020395 

HQ315878 

SA79 

HQ020396 

HQ315879 

SA85 

HQ020397 

HQ315880 

SA98 

HQ020398 

HQ315881 

Colombia,  Amazonas,  Leticia 

C08921 

HQ020382 

- 

(-04.1158,-69.9522) 

C09041 

HQ020383 

HQ315867 

An.  nuneztovari  C 

Colombia,  Antioquia,  Zaragoza 

Za-106 

HQ020399 

- 

(07.4847,  -74.8679) 

Za-110 

HQ020400 

- 

Za-130 

HQ020401 

- 

Za-131 

HQ020402 

- 

Za-158 

HQ020403 

- 

Za-207 

HQ020404 

- 

Za-208 

HQ020405 

- 

Colombia,  Norte  de  Santander,  Tibu 

NS-17 

HQ020388 

HQ315872 

(08.6403,  -72.7372) 

NS-62 

HQ020391 

HQ315875 

NS-71 

HQ020392 

HQ315876 

NS-72 

HQ020393 

HQ315877 

NS-130 

HQ020384 

HQ315868 

NS-135 

HQ020385 

HQ315869 

NS-139 

HQ020386 

HQ315870 

NS-144 

HQ020387 

HQ315871 

NS -243 

HQ020389 

HQ315873 

NS-319 

HQ020390 

HQ315874 

Outgroups 

Anopheles  oswaldoi  B 

Colombia,  Antioquia,  Nechi 
(08.1100,  -74.7670) 

C09000 

- 

HQ315866 

Anopheles  rangeli 

Brazil,  Mato  Grosso 
(-10.2256,  -54.9861) 

BR704-40 

- 

HQ315865 

Anopheles  trinkae 

Ecuador,  Napo,  Sardina  Yacu 

SY205 

- 

HQ315882 

(-00.1666,  -77.0833) 

SY227 

- 

HQ315883 

from  these  two  specimens  and  a  further  17  individuals 
of  An.  nuneztovari  from  the  departments  of  Zaragoza  (n 
=  7)  and  Norte  de  Santander  (n  =  10)  (Table).  Given  the 
morphological  similarities  between^/?,  nuneztovari  and 

An.  dunhami  and  the  reported  presence  of  the  latter  spe¬ 
cies  in  the  Brazilian  Amazon,  DNA  was  also  extracted 
from  abdomens  of  five  archive  specimens  of  An.  dunha¬ 
mi  collected  at  the  type  locality  (Tefe,  Brazil).  Voucher 
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specimens  are  stored  in  the  mosquito  collections  of  the 
Smithsonian  Institution  of  National  Museum  of  Natural 
History  (NMNH),  Washington  DC,  USA  (Table). 

Molecular  methods  -  DNA  template  was  acquired  us¬ 
ing  the  commercially  available  DNeasy  Blood  &  Tissue 
Kit  and  BioSprint  96  DNA  (both  QIAgen®,  Maryland, 
USA).  Amplification  of  rDNA  ITS2  polymerase  chain 
reaction  (PCR)  was  carried  out  using  the  primers  of  Col¬ 
lins  and  Paskewitz  (1996)  and  the  protocol  described  by 
Linton  et  al.  (2001). 

Amplification  of  the  710  bp  fragment  of  the  barcoding 
region  of  the  mitochondrial  COI  gene  was  achieved  us¬ 
ing  the  primer  pair  LCO  1490  and  HCO  2198  developed 
by  Folmer  et  al.  (1994),  using  a  PCR  protocol  optimized 
by  the  Laboratories  of  Analytical  Biology,  Museum  Sup¬ 
port  Center,  Smithsonian  lnstitution/NMNH.  Each  PCR 
contained  1  x  NH4  buffer,  0.5  mM  each  dNTP,  2  mM  of 
MgCl,  0.2  U  of  Taq  polymerase  (BioLine,  Taunton,  MA, 
USA),  0.3  pM  each  primer  and  1  pL  of  DNA  template 
and  was  made  up  to  a  total  volume  of  10  pL  using  ddH,0. 
The  PCR  thermocycler  parameters  included:  a  single 
cycle  at  95°C  for  5  min,  followed  by  34  cycles  of  95°C 
for  30  s,  48°C  for  30  min  and  72°C  for  45  s,  respectively, 
terminating  with  a  72°C  for  5  min  extension  step  and 
a  10°C  hold.  PCR  products  (both  1TS2  and  COI)  were 
visualized  on  1%  agarose  gels,  containing  0.5  mg/mL  of 
ethidium  bromide  and  purified  using  ExoSAP-IT®  (USB 
Corporation,  Cleveland,  Ohio,  USA). 

Sequencing  reactions  were  carried  out  in  both  di¬ 
rections  using  the  Big  Dye  Terminator  Kit®  (PE  Ap¬ 
plied  BioSystems,  Warrington,  England)  on  an  AB1 
3730  automated  sequencer  (PE  Applied  BioSys¬ 
tems).  Sequences  were  edited  using  Sequencher™ 
v.4.8  (Genes  Codes  Corporation,  Ann  Arbor, 
MI)  and  either  aligned  automatically  in  ClustalX 
(Thompson  et  al.  1997)  or  manually  using  MacClade 
v.4.06  (Maddison  &  Maddison  2003).  Sequence  similar¬ 
ities  were  compared  with  those  available  in  GenBank  us¬ 
ing  Basic  Local  Alignment  Search  Tool  (blast. ncbi.nlm. 
nih.gov/Blast.cgi).  Sequence  divergence  was  calculated 
in  MEGA  version  4  (Kumar  et  al.  2008)  using  the  Kimu- 
ra  2-parameter  distance  model  (K2P)  (Kimura  1980).  A 
bootstrap  neighbor  joining  (NJ)  tree  (Saitou  &  Nei  1987) 
was  generated  in  PAUP  version  4.0  (Swofford  2003),  us¬ 
ing  1,000  replicates  of  the  K2P  distance  matrices. 

COI  sequences  generated  formed  part  of  the  ac¬ 
tivities  of  the  Mosquito  Barcoding  Initiative  (project 
leader  YML)  and  both  full  specimen  and  COI  sequence 
data,  including  electropherograms,  are  freely  available 
through  the  Barcode  of  Life  database  website  at  boldsys- 
tems.org.  GenBank  accession  for  sequences  generated  in 
this  study  are  listed  in  Table. 

RESULTS 

ITS2  sequences  (n  =  5)  and  COI  DNA  barcodes  (n 
=  4)  were  generated  from  five  topotypic  archive  speci¬ 
mens  of  An.  dunhami  from  Tefe.  ITS2  sequences  from 
19  Colombian  Anopheles,  morphologically  identified 
as  An.  nuneztovari  using  the  keys  of  Faran  (1980)  and 
Suarez  et  al.  (1988)  revealed  two  distinct  ITS2  sequenc¬ 


es,  differing  at  17  bases  [1  transition  (G/A),  1  transver¬ 
sion  (C/G)  and  15  indels]  (Fig.  1). 

The  most  frequent  ITS2  sequence  was  shared  by 
seven  specimens  from  Zaragosa  (Antioquia)  and  10 
specimens  from  Tibu  (Norte  de  Santander),  morpho¬ 
logically  identified  as  An.  nuneztovari  s.l.,  whereas  a 
different  sequence  was  shared  by  two  specimens  col¬ 
lected  in  Kilometro  12,  Leticia  (Amazonas,  Colombia). 
The  most  frequent  sequence  was  identical  to  previously 
published  sequences  of  An.  nuneztovari  C  from  Colom¬ 
bia  (Sierra  et  al.  2004)  (AY028081-124,  AY028128), 
Brazil  (Marrelli  et  al.  2005)  (AF461749)  and  Venezuela 
(Fritz  et  al.  1994)  (L22462). 

The  two  specimens  from  Kilometro  12,  Leticia 
showed  100%  identity  with  topotypic  An.  dunhami  se¬ 
quences  and  three  published  An.  dunhami  sequences 
from  Parintins  (EU848337,  EU848338  and  EU848344) 
(Calado  et  al.  2008).  However,  sequence  differences  were 
noted  between  the  topotypic  and  Colombian  specimens 
and  the  remaining  six  available  An.  dunhami  sequences 
(generated  from  clones)  reported  from  Parintins.  Topo¬ 
typic  and  Colombian  An.  dunhami  shared  99.7%  identity 
with  EU848340  and  EU848341  (1/360  variant  bases), 
99.4%  identity  with  EU848339,  EU848343,  EU84845, 
same  haplotype  (2/360  variant  bases)  and  98.9%  with 
EU848342  (4/360  variant  bases)  (Fig.  1). 

MtDNA  COI  DNA  barcode  sequences  (658  bp  with¬ 
out  primers)  were  generated  for  specimens  previously 
identified  by  ITS2  sequence  data  as  follow:  topotypic  An. 
dunhami  from  Brazil  (n  =  4)  and  Colombia  (Kilometro 
12,  Leticia)  (n  =  1)  and  An.  nuneztovari  C  from  Norte  de 
Santander  (n  =  10).  COI  barcode  data  were  also  gener¬ 
ated  for  An.  trinkae  (n  =  2),  An.  oswaldoi  B  (n  =  1)  and 
An.  rangeli  (n  =  1),  as  these  species  can  be  morphologi¬ 
cally  mistaken  for  An.  dunhami  or  An.  nuneztovari  s.l.  in 
the  female  adult  stage  (Faran  1980). 

The  bootstrap  NJ  tree  (Fig.  2),  based  on  1,000  rep¬ 
licates  of  the  K2P  COI  distance  matrices  (Kimura 
1980),  reaffirmed  the  specific  identity  of  topotypic 
and  Colombian  An.  dunhami  and  its  genetic  differen¬ 
tiation  from  other  closely  related  taxa  including  An. 
nuneztovari  C,  An.  trinkae,  An.  oswaldoi  B  and  An. 
rangeli  obtained  by  ITS2.  Mean  intra-specific  K2P 


12222222222333333333 
292247888869011225555 
492339 01 234084 54 5345 G 

An.  dunhami  (n*5),  Bra.,  Tef6  t - aagagagaaga - 

An.  dunhami  (n*2),  Col.,  Amazonas  . . . . 

An.  dunhami  {EU848337-38, 44) ,  Bra.,  Par.  . . . . 

An.  dunhami  (EU848339,  43,  45) ,  Bra.,  Par.  . . . — - - 

An.  dunhami  (EU848340),  Bra.,  Par.  . - g . . 

An.  dunhami  (EU848341),  Bra.,  Par.  . . g - 

An.  dunhami  (EU848342),  Bra.,  Par.  c . — c... - 

An.  nuneztovari  (n-17),  Col.,  N.S.,  Zar.  .ggg. - .gc.tgagag 

Fig.  1:  alignment  showing  differing  bases  of  second  internal  tran¬ 
scribed  spacer  sequences  of  Colombian  Anopheles  dunhami  (546  bp; 
n  =  2)  and  Anopheles  nuneztovari  C  (549  bp;  n  =  17)  generated  in  this 
study  in  comparison  to  topotypic  An.  dunhami  from  Tefe,  Brazil  (546 
bp;  n  =  5)  and  other  available  sequences  available  in  GenBank  (360 
bp;  EU848337-42)  (Calado  et  al.  2008).  Bra:  Brazil;  Col:  Colombia; 
Par:  Parintins;  N.S.:  Norte  de  Santander;  Zar:  Zaragoza. 
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divergence  values  were  0.007,  0.004  and  0.010  in  An. 
dunhami  (n  =  5),  An.  nuneztovari  C  (n  =  10)  and  An. 
trinkae  (n  =  2),  respectively.  The  mean  K2P  value  be¬ 
tween  An.  dunhami  and  An.  nuneztovari  was  0.025, 
An.  dunhami  and  An.  trinkae  was  0.093,  An.  nunezto¬ 
vari  and  An.  trinkae  was  0.085. 

DISCUSSION 

This  study  provides  DNA  sequence  data  for  topotypic 
specimens  of  An.  dunhami  for  the  first  time,  thus  estab¬ 
lishing  the  genetic  identity  of  the  species.  Correlation  of 
the  ITS2  rDNA  and  COI  barcodes  of  Colombian  speci¬ 
mens  from  the  Colombian  Amazon  with  topotypic  An. 
dunhami  confirms  the  first  record  of  An.  dunhami  outside 
Brazil  and  comprises  a  new  record  for  Colombia.  That 
An.  dunhami  can  be  morphologically  misidentified  with 
vector  species,  including  An.  nuneztovari  s.l.  (as  in  this 
study),  An.  trinkae  and  An.  oswaldoi  s.l.,  is  of  significance 
to  vector  control  personnel  in  Colombia,  as  it  could  result 
in  wasted  resources,  with  unnecessary  control  measures 
mistakenly  being  applied  against  this  non-vector  species. 

Topotypic  ITS2  sequences  of  An.  dunhami  were  less 
than  1%  dissimilar  to  those  reported  by  Calado  et  al. 
(2008).  The  two  remaining  GenBank  entries  labelled  as 
An.  dunhami  are  problematic:  U9236  (Tefe)  (Lounibos 
et  al.  1998)  and  AF462378  (Acre,  Brazil)  (Marrelli  et  al. 
2005).  As  highlighted  by  Calado  et  al.  (2008),  signifi- 


I  An.  trinkae 


C09000osvv6  An.  oswaldoi  6 
BR704_40ran  An.  rangeli 
SY227trink 
$Y205trink 
NS_319nun 
NS_135nun 
NS  130nun 
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Fig.  2:  bootstrap  neighbor-joining  tree  constructed  using  1,000  repli¬ 
cates  of  Kimura  two-parameter  distance  model  genetic  distance  matri¬ 
ces  of  cytochrome  c  oxidase  unique  haplotypes  (680  bp)  belonging  to  23 
specimens  (only  unique  haplotype  are  shown)  of  Anopheles  nuneztovari 
[Norte  de  Santander,  Colombia  (n=  10),  GenBankHQ3 1 5868-HQ3 15877], 
Anopheles  dunhami  [Kilometro  12,  Amazonas,  Colombia  (n  =  1), 
HQ315867;  Tefe,  Amazonas,  Brazil  (n  =  4),  HQ315878-HQ315881; 
EU848330;  EU848330  (same  haplotype  as  EU848329),  EU848331, 
EU848332  (Calado  et  al.  (2008)].  Out-groups:  Anopheles  trinkae  (Sar- 
dina  Yacu,  Napo  State,  Ecuador,  HQ315882-HQ3 15883),  Anopheles 
oswaldoi  B  (HQ315866)  wad  Anopheles  rangeli  (HQ315865). 


cant  differences  exist  between  the  GenBank  submission 
(U9236)  and  the  published  sequences  in  the  paper  of 
Lounibos  et  al.  (1998)  and  the  submitted  sequence  shows 
15  base  differences  (all  indels)  from  our  An.  dunhami 
sequences.  An  author  of  the  present  study  (JEC)  and  co¬ 
author  of  Lounibos  et  al.  (1998)  confirms  that  manual  se¬ 
quencing  errors  are  most  likely  to  account  for  these  dis¬ 
crepancies.  GenBank  entry  AF462378  of  Marrelli  et  al. 
(2005)  was  more  similar  (96%)  to  An.  oswaldoi  s.s.  from 
Brazil,  Sao  Paulo  (AF425915),  than  topotypic  An.  dun¬ 
hami  86%  (68/496  variant  bases). 

Although  the  two  An.  dunhami  specimens  detected 
were  field  identified  by  females  as  An.  nuneztovari  s.l., 
these  taxa  are  easily  differentiated  using  either  ITS2  or  COI 
sequence  data  or  male  genitalia  (Calado  et  al.  2008).  An. 
dunhami  and  An.  nuneztovari  C  share  only  96.9%  identity 
at  ITS2  (Fig.  1)  and  97.5%  at  COI  (Fig.  2).  An.  dunhami 
ITS2  sequences  generated  in  this  research  share  78-80% 
identity  with  An.  trinkae  sequences  in  GenBank  (U92346, 
U92355  and  Y09075)  (Fritz  1998,  Lounibos  et  al.  1998). 

Contrary  to  the  low  level  of  intra-specific  COI  se¬ 
quence  divergence,  the  mean  inter-specific  K2P  distance 
between  An.  dunhami.  An.  nuneztovari  and  An.  trinkae 
was  0.067  (range  0.025-0.093).  Comparable  results  were 
reported  in  other  barcode  studies  on  mosquitoes,  i.e.  0.058 
between  Anopheles  fluviatilis  and  Anopheles  minimus 
(Kumar  et  al.  2007)  and  10.4%  (range  0.2-17.2%)  for  nine 
genera  of  Culicidae  (Cywinska  et  al.  2006).  In  contrast 
to  the  ITS2  dataset,  using  mtDNA  markers,  An.  dunha¬ 
mi  and  An.  nuneztovari  appear  most  closely  related  with 
2.5%  sequence  divergence,  but  significantly  higher  than 
intra-specific  comparisons. 

Herein  we  verified  the  presence  of  An.  dunhami  in  Co¬ 
lombia  for  the  first  time  using  both  ITS2  sequences  and 
COI  barcodes.  Although  An.  dunhami  has  not  been  impli¬ 
cated  in  malaria  transmission  (Peyton  1993,  Lounibos  et 
al.  1998),  adult  females,  the  stage  most  readily  collected  in 
epidemiological  studies,  can  be  easily  misidentified  as  An. 
nuneztovari  s\.,An.  oswaldoi  s.l.  or  An.  trinkae,  all  of  which 
are  proven  natural  vectors  of  malaria  in  South  America 
(Faran  1980,  Hayes  et  al.  1987,  Herrera  et  al.  1987,  Olano 
et  al.  2001,  Quinones  et  al.  2006,  Gutierrez  et  al.  2009).  It  is 
possible  that  An.  dunhami  may  also  be  a  competent  vector; 
however,  its  role  in  local  malaria  transmission  is  masked 
by  the  presence  of  other  morphologically  similar  species, 
which  are  traditionally  known  as  malaria  vectors.  Knowl¬ 
edge  of  mosquito  species  present  in  vector  control  regions  is 
important.  In  Colombia,  local  vector  control  strategies  vary 
according  to  the  vector  species  present  and  this  in  turn  is 
driven  by  the  entomological  surveys,  thus  it  is  now  impera¬ 
tive  that  morphological  keys  in  Colombia  are  updated  to 
include  An.  dunhami  and  to  determine  novel  morphological 
characters  to  differentiate  this  species  from  An.  nuneztovari 
in  the  adult  female  stage  if  possible.  It  seems  likely  that  this 
species  will  also  be  present  in  the  Amazonian  Region  of 
Peru  and  that  its  eventual  known  distribution  could  include 
a  large  portion  of  the  Amazon  Basin.  Topotypic  ITS2  and 
COI  sequences  of  An.  dunhami  from  Tefe  generated  in  this 
study  will  facilitate  the  identification  of  this  species  across 
its  geographical  range. 
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